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INTRODUCTION 


This  study  is  in  response  to  a request  of  the  Commander,  Scranton 
Army  Ammunition  Plant,  to  determine  if  a proposed  replacement  yard  crane 
system  can  satisfactorily  perform  the  material  handling  task.  The  over- 
head rails  of  the  existing  crane  are  becoming  structurally  inadequate. 
Upon  studying  the  physical  layout  of  the  yard  and  the  structural  integ- 
rity of  the  soil  and  supporting  columns  of  the  present  system,  the  Army 
Corps  of  Engineers  has  recommended  tearing  out  the  present  system  and 
replacing  it  with  a much  slower  traveling  gantry  crane.  Since  a sub- 
stantial amount  of  money  is  involved,  the  Commanding  Officer  at  the 
Scranton  Plant  has  requested  the  development  of  a prove-out  computer- 
ized simulation  model.  In  response  to  this  request,  the  Logistics 
Division  of  the  Systems  Analysis  Directorate  developed  a FORTRAN  based 
computerized  model  described  in  the  balance  of  this  note. 

PURPOSE 

a.  Develop  and  provide  an  analytical  model  which  can  be  used  to 
evaluate  the  capability  of  proposed  alternative  crane  systems  to  per- 
form the  material  handling  task  at  the  Scranton  Army  Ammunition  Plant. 

b.  Determine  if  the  new  yard  crane  proposed  by  the  Army  Corps  of 
Engineers  can  satisfactorily  perform  the  material  handling  task  at  the 
Scranton  Army  Ammunition  Plant. 

METHODOLOGY 

The  FORTRAN  program  developed  for  this  simulation  model  is  event- 
oriented  with  the  capability  to  take  samples  from  the  model  at  discrete 
uniform  time  intervals.  There  are  two  classes  of  events  in  this  model: 

(1)  independent  events,  those  events  the  crane  has  no  control  over,  and 

(2)  dependent  events,  those  events  the  crane  creates  in  responding  to 
the  independent  events.  The  independent  events  are  as  follows: 


1. 

A charging  call  from  Feeder 

or  Line 

No. 

2. 

A charging  call  from  Feeder 

or  Line 

No. 

3. 

A charging  call  from  Feeder 

or  Line 

No. 

4. 

The  arrival  of  a 5-1/4  inch 

heat. 

5. 

The  arrival  of  a 6 inch  heat. 

6. 

The  arrival  of  a 7-3/8  inch 

heat. 

7. 

The  morning  coffee  break. 

8. 

The  lunch  break. 
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9.  The  afternoon  coffee  break 


10.  The  shift  break. 

Specific  distributions  have  been  derived  for  the  length  of  time  between 
repeated  occurrences  of  the  same  event  (see  APPENDIX  B) . 

The  dependent  events  of  preparing  and  loading  charges  into  the 
feeders  and  unloading  and  stacking  billets  in  open  storage  bays  consisted 
of  7 basic  tasks  as  follows: 

1.  Picking  billets  up. 

2.  Setting  billets  down. 

3.  Squaring  billets  up. 

4.  Swinging  billets  into  position. 

5.  Breaking  for  coffee,  lunch  and  shift  turnover. 

6 . Traveling . 

7.  Idleness. 

Specific  time  distributions  for  each  of  the  first  5 tasks  were  derived 
for  each  of  the  unique  operations  in  which  these  tasks  occurred.  The 
standard  Newtonian  Equations  of  Motion  were  employed  to  model  the  crane’s 
travel.  These  equations  were  further  supplied  with  data  generated  from 
acceleration,  deceleration,  and  top  velocity  distributions.  These  dis- 
tributions allowed  for  operator-to-operator  speed  variation.  Break  times 
were  modeled  as  distributions;  however,  the  management  cautioned  that 
these  periods  could  be  eliminated  if  a high  level  of  impending  activity 
warranted  it.  Idleness  was  the  residual.  If  the  crane  had  nothing  to 
do,  it  was  idle.  (See  APPENDIX  B for  a review  of  the  previously  described 
inputs.)  These  7 basic  crane  tasks  were  monitored  during  the  simulation 
to  enable  the  crane's  consumption  of  time  to  be  reviewed. 

The  basic  structure  of  this  model  consists  of  a 1 dimensional 
array  which  stores  the  next  time  each  of  the  10  independent  events  will 
occur,  and  a clock  which  is  event-incremented.  These  events  include  all 
the  independent  events  and  the  various  tasks  of  the  dependent  events. 

Each  time  the  clock  is  incremented  after  completing  a dependent  event, 
the  next  independent  event  is  checked  to  see  if  it  has  occurred.  If  an 
independent  event  has  occurred,  it  is  processed  and  a new  time  is  generated 
to  determine  when  this  event  will  next  occur.  If  an  independent  event 
has  not  occurred,  the  program  continues  processing  the  crane's  dependent 
events.  If  there  are  not  any  dependent  events  to  process,  the  model 
advances  to  processing  the  next  independent  event. 
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ASSUMPTIONS 


When  constructing  a mathematical  model,  some  assumptions  or  ground 
rules  must  be  formed  to  define  a base  for  the  study.  Further,  some 
assumptions  must  be  made  on  noncritical  entities  to  avoid  creating  more 
structure  than  is  necessary  to  create  a desired  or  adequate  level  of 
realism.  A high  level  of  resolution  of  noncritical  activities  cannot  be 
economically  justified.  These  assumptions  (except  the  first  one)  can 
easily  be  altered  within  the  computerized  model.  They  are  as  follows: 

1.  Only  one  crane  will  be  in  operation  at  a time. 

2.  100  percent  reliability  is  required  (at  least  one  crane  will  be 
operable  at  all  times). 

3.  The  billets  arrive  in  groups  of  cars  carrying  one  heat  rather 
than  in  individual  cars  strung  out  over  an  extended  period  of  time,  thus 
creating  an  intermixing  of  heats. 

4.  The  heats  will  be  unloaded  in  car  segments  with  the  cars  being 
individually  positioned  with  an  offset  from  the  center  line  of  the  receiv- 
ing bay. 


5.  Billets  are  stored  and  charged  or  loaded  into  the  feeders  in 
heats  which  are  assumed  to  be  groups  of  about  175  billets  each. 

6.  Bays  are  assigned  usage  priorities  and  a maximum  load  capacity 
for  the  given  type  of  billet  (5-1/4,  6,  or  7-3/8  inches)  which  the  bay 
can  handle. 

7.  At  least  one  charge  for  all  three  feeders  will  be  ready  at  all 
times  and  will  be  residing  in  a work  area  where  the  billets  are  assimilated 
into  squared  groupings  ready  for  loading  or  charging  into  the  feeder. 

8.  When  answering  a feeder  request  for  a charge,  the  crane,  as  soon 
as  it  releases  any  billets  it  might  be  carrying,  will  load  the  feed  table. 
Another  charge  will  be  immediately  prepared  for  that  feeder  unless  another 
feeder  calls  for  a charge;  in  which  case  the  crane  will  load  that  feeder 
and  then  proceed  to  prepare  charge (s)  for  any  feeders  not  having  a charge 
in  reserve. 

9.  Breaks  occur  throughout  the  plant.  Therefore,  no  crane  activity 
occurs  during  any  of  the  breaks. 

10.  Crane  travel  is  uninterrupted.  Lift  trucks,  etc.  don’t  get  in 
the  way  of  the  crane. 

11.  The  assumed  level  of  manufacturing  activity  will  be  for  mobil- 
ization. 
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12.  The  crane  operator  has  sufficient  skill  and  is  unencumbered 
by  union  control  or  management  policy  to  allow  him  to  operate  the  crane 
at  maximum  speed.  He  could,  thus,  start  the  X and  Y movements*,  simul- 
taneously, before  completely  clearing  a railroad  car  or  bay  stock-pile 
containing  a load  of  billets,  i.e.,  not  be  confined  to  slow  rectangular 
motion. 


13.  When  computing  travel  times  and  distances,  the  X and  Y center 
line  coordinates  of  the  bays,  work  areas,  feeders,  and  railroad  cars 
are  used  to  calculate  the  distance  between  these  objects. 

DISCUSSION 

After  making  some  preliminary  test  runs  to  validate  the  model,  two 
runs  were  made  to  verify  the  adequacy  of  the  200-feet-per -minute  top 
velocity  of  the  proposed  new  crane.  Both  runs  used  the  previously  des- 
cribed priority  and  bay  layout  in  APPENDIX  A plus  the  numeric  inputs  as 
described  in  APPENDIX  B. 

To  review  the  major  numeric  items,  the  top  velocity  in  the  X and  Y 
direction  was  entered  as  a uniform  distribution  ranging  between  150  and 
200  feet-per-minute.  This  range  allowed  for  operator-to-operator  and 
circumstance- to-circumstance  variation.  Additionally,  the  acceleration 
and  deceleration  rates  in  the  X and  Y directions  were  entered  as  a con- 
stant of  1-f oot-per-second  squared.  Since  200  feet-per-minute  is  equiva- 
lent to  only  2.23  miles  per  hour,  it  was  generally  believed  that  accelera- 
tion and  deceleration  to  and  from  this  level  would  take  place  at  the  maximum 
rate  of  1-f oot-per-second  squared.  Further,  at  this  maximum  acceleration 
rate,  the  small  distance  of  5.56  feet  is  required  to  achieve  this  top  velocity 
of  200  feet-per-minute  when  starting  from  a standstill.  These  runs  were 
simulated  for  a period  in  excess  of  1 year  under  steady  state  conditions. 
Additionally,  these  runs  did  not  have  any  unload  queue  build-up  and  all 
the  bays  in  the  storage  yard  had  an  inventory  of  175  billets — 1 heat. 

Run  1 represented  a fairly  good  way  of  operating.  While  not  optimal, 
the  bay  layout  and  priority  scheme  used  was  believed  to  be  near  optimal 
if  management  kept  the  work  or  staging  areas  in  their  present  location 
as  exhibited  on  the  layout  in  APPENDIX  A.  Under  these  conditions,  the 
crane  was  able  to  handle  billets  at  a fast  enough  rate  to  enable  it  to 
be  idle  27%  of  the  average  day  with  an  additional  break  time  (2  coffee 
and  1 lunch  and  1 shift)  idleness  of  nearly  18%.  Additionally,  the  feed 
tables  had  to  wait  an  average  of  only  2-3/4  minutes  for  a recharge.  Once, 
feeder  table  number  2 had  to  wait  9.19  minutes  for  a recharge;  however, 
the  histogram  plots  of  wait  times  reveal  that  this  situation  has  very 


The  X direction  is  the  direction  of  travel  parallel  to  the  crane  rails.  The 
Y direction  is  perpendicular  to  the  rails.  See  the  bay  layout  in  APPENDIX  A. 
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little  probability  of  occurrence.  These  histograms  showed  that,  in  gen- 
eral, the  feed  table  recharge  wait  time  was  5 minutes  or  less  95%  of  the 
time,  and  about  6-1/4  minutes  or  less  99%  of  the  time.  It  was  further 
observed  in  some  of  the  validation  runs  that  moving  the  work  area  or 
staging  area  directly  across  the  railroad  tracks  to  the  open  space  pre- 
sently on  the  layout  and  utilizing  the  present  work  area  for  storage 
space  could  significantly  reduce  the  feed  table  recharge  wait  time. 

Input  data  for  run  2 was  the  same  as  for  run  1 except  the  bay 
priorities  were  inverted,  i.e.,  the  bays  farther  away  from  the  feed  tables 
were  used  the  most.  This  setup  represents  a very  poor  operating  pro- 
cedure. Under  these  conditions,  the  crane  was  idle  only  about  7%  of  the 
typical  day.  The  average  feed  table  wait  time  increased  about  1-3/4 
minutes  to  yield  an  average  of  4-1/2  minutes.  Additionally,  the  maxi- 
mum wait  times  observed  increased  along  with  an  exhibited  downward  shift 
in  the  wait  time  histograms.  However,  under  these  adverse  operating 
conditions,  the  crane  still  performed  adequately. 

Finally,  a third  run  was  made  to  observe  how  the  crane  would  respond 
to  a shock  load.  Run  3 was  initiated  at  the  end  of  run  1 using  the 
ending  values  of  the  simulation  variables  obtained  at  the  end  of  run  1 
with  the  exception  of  an  addition  of  1 month’s  consumption  of  inventory 
being  placed  In  the  unload  queues.  It  took  the  crane  770  hours  or  about 
1-1/4  months  to  empty  the  unload  queues  and  thus  return  to  steady  state 
conditions.  While  the  crane  was  working  off  this  extra  load,  business 
was  being  conducted  as  usual.  It  would  thus  appear  that  this  response 
is  more  than  adequate  since  management  does  not  expect  any  shocks  of 
this  magnitude  to  be  imposed  on  the  system. 

CONCLUSION 

The  computer  runs  reviewed  In  the  DISCUSSION  section  amply  illustrate 
that  the  200-f eet-per-rainute  top  velocity  gantry  crane  proposed  by  the 
Army  Corps  of  Engineers  should  more  than  adequately  be  able  to  handle 
the  proposed  mobilization  workload.  Further,  it  is  apparent  that  the 
priority  arrangement  of  bays  and  the  location  of  the  work  or  staging 
areas  can  critically  affect  the  crane  workload.  It  is  strongly  urged 
that  this  model  be  exercised  to  optimize  the  preceding  2 parameters . 

Finally,  the  various  runs  made  for  this  study  along  with  a critical 
review  by  the  Scranton  management  indicates  that  the  model  is  very  realistic. 
Thus,  the  objectives  of  this  study  as  listed  in  the  PURPOSE  section  of 
this  report  have  been  satisfactorily  completed. 


Next  page  is  blank. 
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APPENDIX  B 


CARD  LAYOUTS  AND  DATA 


A.  CONTROL  CARDS 

First  Card 

Col  1 (Format  II)  Enter  a zero  in  this  column  or  leave  blank  if 
the  event  and  inventory  listing  is  desired  for  the  length  of  time  con- 
strained by  the  times  entered  in  fields  2 (cols  7-12)  and  3 (cols  13-18) 
of  this  card.  Otherwise,  enter  a 1 in  this  column  to  suppress  the 
listing. 

Col  2 (Format  II)  Enter  a zero  in  this  column  or  leave  blank  if 
the  bay  trace  listing  is  desired  for  the  length  of  time  constrained  by 
the  times  entered  in  fields  2 (cols  7-12)  and  3 (cols  13-18)  of  this 

card.  Otherwise,  enter  a 1 in  this  column  to  suppress  the  listing. 

Col  3 (Format  II)  Enter  a zero  in  this  column  or  leave  blank  if 

the  time  listing  is  desired  for  the  length  of  time  constrained  by  the 

times  entered  in  fields  2 (cols  7-12)  and  3 (cols  13-18)  of  this  card. 
Otherwise,  enter  a 1 in  this  column  to  suppress  the  listing. 

Col  4 (Format  II)  Enter  a zero  in  this  column  or  leave  blank  if  the 
distance  listing  is  desired  for  the  length  of  time  constrained  by  the 
times  entered  in  fields  2 (cols  7-12)  and  3 (cols  13-18)  of  this  card. 
Otherwise,  enter  a 1 in  this  column  to  suppress  the  listing. 

Col  5 (Format  II)  Enter  a zero  in  this  column  or  leave  blank  if  a 
50  observation  listing  of  generations  of  the  stochastic  data  entered  in 
this  program  is  desired.  Otherwise,  enter  a 1 in  this  column  to  suppress 
the  listing. 

Col  6 (Format  II)  Enter  a zero  in  this  column  or  leave  blank  if  a 
time  vs.  inventory  plot  of  the  bays  is  desired  and  a histogram  plot  of 
the  work  areas  is  desired.  Otherwise,  enter  a 1 in  this  column  to  suppress 
this  listing. 

Cols  7-12  (Format  F6.0)  Enter  (in  minutes)  the  time  when  it  is  desired 
to  start  listing  the  entries  of  columns  1-4  above. 

Cols  13-18  (Format  F6.0)  Enter  (in  minutes)  the  time  when  it  is 
desired  to  stop  listing  the  data  in  the  preceding  defined  lists. 

Cols  19-24  (Format  F6.0)  Enter  (in  minutes)  the  time  when  it  is  desired 
to  stop  the  simulation. 

Cols  25-26  (Format  F2.0)  Enter  (in  minutes)  the  length  of  time 
desired  between  intervals  when  taking  observations  of  the  simulation. 
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Cols  27-30  (Format  F4.1)  Enter  (in  feet)  according  to  the  origin 
defined  in  the  yard  layout,  the  X coordinate  of  the  center  line  of  the 
crane  magnet  when  starting  the  simulation. 

Cols  31-34  (Format  F4.1)  Enter  (in  feet)  the  Y coordinate  as 
defined  in  the  preceding  field. 

Cols  35-44  (Format  110)  Enter  the  value  initially  assigned  to  the 
seed  of  the  uniform  (0  to  1)  random  number  generator.  The  last  value 
of  the  seed  will  be  punched  out  at  the  end  of  the  run.  This  seed  may 
then  be  used  for  a follow  on  run,  thus  enabling  the  user  to  make  a 
replication  over  a different  sequence  of  random  numbers.  (A  suggested 
initial  seed  is  435459  - IBM's  RANDU) . 

Cols  45-46  (Format  12)  Enter  the  type  of  material  desired  to  be 
unloaded  first  - 0 (zero)  or  blank  = no  particular  choice,  1 = 5-1/4, 
2=6,  3 = 7-3/8.  If  something  other  than  0 (zero)  or  blank  is  entered 
in  this  field,  some  of  the  indicated  type  of  billets  must  reside  in 
the  unload  heat  queue. 

Cols  47-48  (Format  12)  Enter  the  bay  number  where  the  billets  are 
to  be  unloaded.  Leave  blank  or  enter  a 0 (zero)  if  it  is  desired  to  let 
the  program  select  the  bay. 

Cols  49-50  (Format  12)  Enter  the  heat  level  of  the  preceding  bay 
which  is  receiving  the  billets  being  unloaded.  Leave  this  field  blank 
or  enter  a 0 (zero)  if  it  is  desired  to  let  the  program  select  the  heat 
level . 

Cols  51-52  (Format  12)  Enter  the  number  of  the  bay  currently  supply- 
ing the  5-1/4  inch  feed  table  with  billets.  Leave  blank  or  enter  a 
0 (zero)  if  it  is  desired  to  let  the  program  perform  this  task. 

Cols  53-54  (Format  12)  Enter  the  number  of  the  bay  currently  supply- 
ing the  6 inch  feed  table  with  billets.  Leave  blank  or  enter  a 0 (zero) 
if  it  is  desired  to  let  the  program  perform  this  task. 

Cols  55-56  (Format  12)  Enter  the  number  of  the  bay  currently  supply- 
ing the  7-3/8  inch  feed  table  with  billets.  Leave  blank  or  enter  a 
0 (zero)  if  it  is  desired  to  let  the  program  perform  this  task. 

Cols  57-58  (Format  12)  Enter  the  heat  level  of  the  bay  currently 
providing  the  5-1/4  inch  billets.  Leave  this  field  blank  or  enter  a 
0 (zero)  if  it  is  desired  to  let  the  program  select  this  heat  level. 

Cols  59-60  (Format  12)  Enter  the  heat  level  of  the  bay  currently 
providing  the  6 inch  billets.  Leave  this  field  blank  or  enter  a 
0 (zero)  if  it  is  desired  to  let  the  program  select  this  heat  level. 
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Cols  61-62  (Format  12)  Enter  the  heat  level  of  the  bay  currently 
providing  the  7—3/8  inch  billets.  Leave  this  field  blank  or  enter  a 
0 (zero)  if  it  is  desired  to  let  the  program  select  this  heat  level. 

Cols  63-64  (Format  12)  Enter  the  minimum  number  of  5-1/4  inch 
billets  required  in  the  squaring-up  charges  work  area. 

Cols  65-66  (Format  12)  Enter  the  minimum  number  of  6 inch  billets 
required  in  the  squaring  up  charges  work  area. 

Cols  67-68  (Format  12)  Enter  the  minimum  number  of  7-3/8  inch 
billets  required  in  the  squaring-up  charges  work  area. 

Cols  69-70  (Format  12)  Enter  the  number  of  5-1/4  inch  billets  in 
a standard  charge  or  loading  of  the  feeder. 

Cols  71-72  (Format  12)  Enter  the  number  of  6 inch  billets  in  a 
standard  charge  or  loading  of  the  feeder. 

Cols  73-74  (Format  12)  Enter  the  number  of  7-3/8  inch  billets 
in  a standard  charge  or  loading  of  the  feeder. 

Cols  75-76  (Format  12)  Enter  the  maximum  number  of  5-1/4  inch 
billets  the  crane  can  carry  to  the  feeder. 

Cols  77-78  (Format  12)  Enter  the  maximum  number  of  6 inch  billets 
the  crane  can  carry  to  the  feeder. 

Cols  79-80  (Format  12)  Enter  the  maximum  number  of  7-3/8  inch 
billets  the  crane  can  carry  to  the  feeder. 

Second  Card 


Cols  1-6  (Format  F6.0)  Enter  (in  feet),  according  to  the  origin 
defined  in  the  yard  layout,,  the  X coordinate  of  the  center  line  of  the 
5-1/4  inch  billet  feeder. 

Cols  7-12  (Format  F6.0)  Enter  (in  feet)  the  Y coordinate  as  defined 
in  the  preceding  field. 

Cols  13-18  (Format  F6.0)  Enter  (in  feet),  according  to  the  origin 
defined  in  the  yard  layout,  the  X coordinate  of  the  center  line  of  the 
6 inch  billet  feeder. 

Cols  19-24  (Format  F6.0)  Enter  (in  feet)  the  Y coordinate  as 
defined  in  the  preceding  field. 

Cols  25-30  (Format  F6.0)  Enter  (in  feet),  according  to  the  origin 
defined  in  the  yard  layout,  the  X coordinate  of  the  center  line  of  the 
7-3/8  inch  billet  feeder. 
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Cols  31-36  (Format  F6.0)  Enter  (in  feet)  the  Y coordinate  as  defined 
in  the  preceding  field. 

Cols  37-42  (Format  F6.0)  Enter  (in  feet),  according  to  the  origin 
defined  in  the  yard  layout,  the  X coordinate  of  the  center  line  of  the 
5-1/4  inch  billet  work  area  used  for  squaring-up  a group  of  billets  to 
be  loaded  into  the  5-1/4  inch  billet  feeder. 

Cols  43-48  (Format  F6.0)  Enter  (in  feet)  the  Y coordinate  as  defined 
in  the  preceding  field. 

Cols  49-54  (Format  F6.0)  Enter  (in  feet),  according  to  the  origin 
defined  in  the  yard  layout,  the  X coordinate  of  the  center  line  of  the 
6 inch  billet  work  area  used  for  squaring-up  a group  of  billets  to  be 
loaded  into  the  6 inch  billet  feeder. 

\ 

Cols  55-60  (Format  F6.0)  Enter  (in  feet)  the  Y coordinate  as  defined 
in  the  preceding  field. 

Cols  61-66  (Format  F6.0)  Enter  (in  feet),  according  to  the  origin 
defined  in  the  yard  layout,  the  X coordinate  of  the  center  line  of  the 
7-3/8  inch  billet  work  area  used  for  squaring-up  a group  of  billets  to 
be  loaded  into  the  7-3/8  inch  billet  feeder. 

Cols  67-72  (Format  F6.0)  Enter  (in  feet)  the  Y coordinate  as  defined 
in  the  preceding  field. 

Cols  73-80  Blank. 

Third  Card 


Cols  1-4  (Format  14)  Enter,  the  number  of  billets  in  the  oldest 
5-1/4  inch  billet  heat  now  residing  in  the  unload  queue. 

Continue  the  preceding  using  fields  of  4 consecutive  columns  until 
the  list  of  5-1/4  heats  queued-up  is  completely  exhausted  or  until  the 
21  heat  is  required.  Since  only  20  heats  are  currently  allowed,  entering 
more  heats  will  require  expanding  the  check  variable  maxque  and  all  the 
arrays  mnemonically  dimensioned  in  terms  of  this  variable  as  shown  in 
subroutine  load. 

Fourth  Card 

Fallow  the  same  procedure  as  outlined  for  the  third  card  except  enter 
data  for  the  6 inch  billets. 

Fifth  Card 

Follow  the  same  procedure  as  outlined  for  the  third  card  except 
enter  data  for  the  7-3/8  inch  billets. 
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Sixth  Card 


Cols  1-8  (Format  F8.0)  Enter  the  future  time  when  independent 
event  number  1 will  next  occur . 

Continue  the  preceding  using  fields  of  8 consecutive  columns  until 
future  event  times  for  all  ten  independent  events  have  been  entered. 

B.  BAY  CARDS 

(One  card  is  required  for  each  bay.  The  last  card  of  this  series  must 
have  -9  punched  in  cols  1-2  and  the  rest  of  the  card  must  be  left  blank.) 

Cols  1-2  (Format  12)  Enter  the  bay  number  (the  bays  must  be  sequen- 
tially numbered) . 

Cols  3-4  (Format  12)  If  this  bay  is  dedicated  to  storing  5-1/4  inch 
billets,  enter  a bay  priority  number  (lower  numbers  are  of  higher  priority) 
Otherwise,  leave  this  field  blank. 

Cols  5-6  (Format  12)  If  this  bay  is  dedicated  to  storing  6 inch 
billets,  enter  a bay  priority  number  (lower  numbers  are  of  higher  priority) 
Otherwise,  leave  this  field  blank. 

Cols  7-8  (Format  12)  If  this  bay  is  dedicated  to  storing  7-3/8  inch 
billets,  enter  a bay  priority  number  (lower  numbers  are  of  higher  priority) 
Otherwise,  leave  this  field  blank. 

Cols  9-14  (Format  F6.0)  Enter  (in  feet),  according  to  the  origin 
defined  in  the  yard  layout,  the  X coordinate  of  the  center  line  of  this 
bay. 


Cols  15-20  (Format  F6.0)  Enter  (in  feet)  the  Y coordinate  as  defined 
in  the  preceding  field. 

Cols  21-25  (Format  15)  Enter  the  maximum  number  of  billets  allowed 
in  this  bay. 

Cols  26-30  (Format  15)  Enter  the  actual  number  of  billets  currently 
inventoried  in  this  bay  for  the  first  heat  level.  Continue  entering, 
using  fields  of  5 consecutive  columns  until  the  billets  of  all  heats 
residing  in  this  bay  have  been  entered. 

NOTE:  All  the  cards  described  up  to  this  point,  except  for  the  first 

26  columns  of  the  first  card  (card  one  of  the  control  cards),  are  punched 
out  at  the  end  of  each  run  with  the  end  of  run  data  entered.  However, 
the  future  event  times  (card  six  of  the  control  cards)  have  been  reduced 
by  the  magnitude  of  the  simulation  stop  time  (columns  19-24  of  control 
card  number  1) . This  setup  will  enable  breaking  long  simulation  runs 
into  segments,  thereby,  enabling  examining  pieces  of  the  run  via  use 
of  some  or  all  of  the  various  long  form  event-by-event  printouts  offered 
in  this  program. 
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C.  STOCHASTIC  DATA  CARDS  CARRY  FIX  POINT  (A  CONSTANT)  AND  RANDOM  (STOCHASTIC) 
DATA  INTO  THE  PROGRAM. 

The  types  of  distributions  available  and  the  definitions  of  inputs 
for  the  various  fields  are  as  follows: 

TYPE  OF  FIELD  FIELD  FIELD  FIELD  FIELD  FIELD  FIELD 

DISTRIBUTION NO.  1 NO.  2 NO.  3 NO.  4 NO.  5 NO.  6 NO.  7 


CONSTANT 

1 

CONSTANT 

UNIFORM 

2 

MIN  OBS 

MAX  OBS 

TRIANGULAR 

3 

MIN  OBS 

MAX  OBS 

MOST 

LIKELY  OBS 

NORMAL 

4 

MIN  OBS 

MAX  OBS 

MEAN 

STANDARD 

DEVIATION 

LOGNORMAL 

5 

MIN  OBS 

MAX  OBS 

MEAN 

STANDARD 

DEVIATION 

GAMMA 

6 

MIN  OBS 

MAX  OBS 

MEAN 

STANDARD 

DEVIATION 

WEIBULL* 

7 

MIN  OBS 

MAX  OBS  SCALE 

PARAMETER 

SHAPE 

PARAMETER 

ERLANG 

8 

MIN  OBS 

MAX  OBS 

MEAN 

NO.  OF 

EXPONETIAL 

DEVIATES 

(EXPONETIAL) 

8 

MIN  OBS 

MAX  OBS 

MEAN 

1 

CHI  SQUARE 

9 

MIN  OBS 

MAX  OBS 

NO. 

DEGREES 

FREEDOM 

*THE  MINIMUM  OBSERVATION  IS  THE  LOCATION  PARAMETER. 
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TYPE  OF  field  field  field  field  field  field 

DISTRIBUTION NO.  1 NO.  2 NO.  3 NO.  4 NO.  5 NO.  6 


4.4. 

BETA 

10 

MIN 

OBS 

MAX  OBS  A 

B 

POISSON 

11 

MIN 

OBS 

MAX  OBS  MEAN 

*** 

PASCAL 

12 

MIN 

OBS 

MAX  OBS  P 

K 

VtVfe  * 

(GEOMETRIC) 

12 

MIN 

OBS 

MAX  OBS  P 

K=  I 

**** 

BINOMIAL 

13 

MIN 

OBS 

MAX  OBS  P 

N 

***** 

HYPERGEOMETRIC 

14 

MIN 

OBS 

MAX  OBS  P 

N 

M 

****** 

HISTOGRAM 

15 

A 

B C 

D 

E 

**F (X)  = G(A+B)X  A_1 

(1-X) 

B-l 

A = Greater 
B = Greater 

than  Zero 
than  Zero 

G(A)G(B) 

G = Gamma  Function 

*** 


F (X) 


**** 

F (X) 


*****  , „ 
F (X) 


(K+X-1)PKQ  x 

X = 0,  1, 

( X ) 

Q = 1-P 

(N)pVX 

X = 0,  1, 

(X) 

Q = 1-P 

(NP)  (NQ  ) 

(X  ) (M-X) 

X = 0,  1, 

(N) 

M-X  = 0,  1, 

(M) 

Q = 1-P 

2... 


2.  . .N 


2.  . .N 
2. . .NQ 


^ ^ ^ ^ 

1.  A - Enter  the  value  of  the  left  hand  boundary  of  probability 
Cell  No.  1. 


2.  B - Enter  the  probability  of  realizing  Cell  No.  1. 

3.  C - Enter  the  value  of  the  right  hand  boundary  of  probability 
Cell  No.  1 and  the  left  hand  boundary  of  probability  Cell  No.  2. 

4.  D - Enter  the  probability  of  realizing  Cell  No.  2. 

5.  E - Enter  the  value  of  the  right  hand  boundary  of  probability 
Cell  No.  2 and  the  left  hand  boundary  of  probability  Cell  No,  3. 

6.  F - Enter  the  probability  of  realizing  Cell  No.  3. 


FIELD 
NO.  7 
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7.  In  Field  No.  1 of  the  next  card  enter  the  value  of  the  right 
hand  boundary  of  probability  Cell  No,  3 and  the  left  hand  boundary  of 
probability  Cell  No,  4. 

8.  Continue  entering  the  elements  of  the  histogram  in  these  7 field 
patterns  using  as  many  cards  as  necessary  until  the  entire  histogram  has 
been  completely  loaded. 

9.  The  next  field  after  the  last  field  used  to  load  the  histogram 
data  must  have  a -999.0  entered  in  it  to  mark  the  end  of  the  histogram 
information. 

The  stochastic  data  card’s  field  layouts  are  as  follows: 

Cols  1-3  (Format  13)  Enter  the  key  number  the  program  uses  to  fetch 
an  observation  from  a given  stochastic  input.  The  data  definitions  of 
the  key  numbers  are  as  follows: 

KEY  NO.  DATA  DEFINITION 


1.  The  distribution  of  resupply  calls  from  Feeder  No,  1.  First  it 
is  necessary  to  calculate  the  time  between  charges  given  the  mobilization 
rate  of  production  per  month  and  the  working  hours  per  month.  Since  the 
calls  of  Feeders  2 and  3 require  similar  information,  some  of  the  necessary 
preliminary  calculations  required  for  them  will  also  be  exhibited. 


ROUND 

GUN  MOB  RATE 

SCRAP 

MULT  WT 

PLANT  INPUT 

RAW  STOCK 

(TONS/MONTH) 

No . /Month 

5% 

In  Lbs 

Tons/Month 

5-1/4" 

6"  7-3/8" 

M107 

155MM 

100,000 

5000 

107 

5618 

7865 

M110 

155MM 

40,000 

2000 

107 

2247 

M437 

175MM 

40,000 

2000 

172 

3612 

3612 

M106 

8 Inch 

40,000 

2000 

220 

4620 

6384 

M404 

8 Inch 

15,000 

750 

224 

1764 

Number  of  Working  Hours  Per  Month,  Using  1976  as  a Base  for  Calculations: 


1. 

Jan 

31  - 5 = 

26 

5. 

May  31  - 6 ■ 

25 

9. 

Sept  30  - 5 = 

25 

2. 

Feb 

29  - 6 = 

23 

6. 

Jun  30  - 4 ■ 

26 

10. 

Oct 

31  - 7 = 

24 

3. 

Mar 

31  - 4 = 

27 

7. 

July  31  - 5 ■ 

26 

11. 

Nov 

30  - 5 = 

25 

4. 

Apr 

30  - 4 - 

26 

8. 

Aug  31  - 5 = 

26 

12. 

Dec 

31  - 5 = 

26 

********Sums  to  305  Working  Days 

Year  (1976) 
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Working  Working  Contact  Contact 

305  Da^£  x Yeaj_  _ 25.42  Days  x 24  Hours  = 610  Hours  Breaks  Total  80  Minutes 

Year  12  Months  Month  Working  Month  8 Contact  Hours 
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Thus,  the  STD  DEV  = (0.15)  (134.27)  = 10.28  also  3 STD  DEV  CUT-OFFS  = 103.44  and  165.10 

1.96 


The  Distribution  of  resupply  calls  from  Feeder  No.  3.  Following  the  calculations  for  Feeder  No. 
yields  the  following: 
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Arrival  of  a 6 Inch  Heat. 
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10.  The  time  from  the  shift  break  to  the  next  occurrence  of  this  same 
event  was  entered  as  a constant  of  480  minutes. 

11.  The  time  required  to  pick  billets  out  of  a railroad  car  was  entered 
as  a Triangular  Distribution  having  a minimum  time  of  0.25  minutes,  a 
maximum  time  of  1.5  minutes,  and  a most-likely  time  of  0.75  minutes. 

This  is  the  time  required  between  stopping  X and  Y motion  to  position 
the  magnet  over  the  car  and  starting  X and  Y motion  to  transport  the  load 
of  billets  to  the  storage  bay. 

12.  The  time  required  to  swing  a load  of  billets  into  position  so  this 
load  can  be  stacked  on  the  pile  in  a storage  bay  was  entered  as  Triangular 
Distribution  having  a minimum  time  of  0,25  minutes,  a maximum  time  of 

1.0  minutes,  and  a most-likely  time  of  0.5  minutes.  Since  this  task  is 
essentially  accomplished  by  a man  on  the  pile,  the  crane’s  travel  to  the 
pile  has  been  completed.  Therefore,  this  time  is  defined  as  the  time 
from  the  end  of  travel — to  the  pile — to  the  start  of  setting  the  load  of 
billets  down  on  the  pile. 

13.  The  time  required  to  set  the  load  of  billets  on  the  storage  pile 
begins  at  the  end  of  swing  time  and  continues  until  the  X and  Y travel 
away  from  the  pile  begins.  This  time  includes  pile  shape  up  time  and 
billet  packing  time.  This  time  was  entered  as  a Triangular  Distribution 
having  a minimum  time  of  0.5  minutes,  a maximum  time  of  2.0  minutes,  and 
a most-likely  time  of  0.85  minutes. 

14.  The  X distance  in  feet  between  the  center  line  of  the  railroad  car 
and  the  receiving  storage  bay  was  entered  as  a histogram  as  follows: 

20%  of  the  time  the  distance  will  be  between  0 and  50  feet,  20%  of  the 
time  the  the  distance  will  be  between  50  and  100  feet,  30%  of  the  time 
the  distance  will  be  between  100  and  150  feet,  and  finally  30%  of  the 
time  the  distance  will  be  between  150  and  200  feet. 

15.  The  number  of  5-1/4  inch  billets  picked  out  of  a railroad  car  per 
unit  pick  was  entered  as  a Triangular  Distribution  having  1.0  as  the 
minimum  number  of  billets,  12.0  as  the  maximum  number  of  billets,  and 

8.0  as  the  most-likely  number  of  billets. 

16.  The  number  of  6 inch  billets  picked  out  of  a railroad  car  per  unit 
pick  was  entered  as  a Triangular  Distribution  having  1.0  as  the  minimum 
number  of  billets,  12.0  as  the  maximum  number  of  billets,  and  7.0  as  the 
most-likely  number  of  billets. 

17.  The  number  of  7-3/8  inch  billets  picked  out  of  a railroad  car  per 
unit  pick  was  entered  as  a Triangular  Distribution  having  1.0  as  the 
minimum  number  of  billets,  8.0  as  the  maximum  number  of  billets,  and 

6.0  as  the  most-likely  number  of  billets. 

18.  The  length  of  the  morning  coffee  break  was  entered  as  a constant  of 

15.0  minutes. 
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19.  The  length  of  the  lunch  break  was  entered  as  a constant  of  35,0 
minutes . 

20.  The  length  of  the  afternoon  coffee  break  was  entered  as  a constant 
of  15.0  minutes. 


21.  The  lepgth  of  the  shift  break  was  entered  as  a constant  of  15.0 
minutes. 

22.  The  time  required  to  pick  billets  off  of  the  storage  pile  in  some 
bay  was  entered  as  a Triangular  Distribution,  having  a minimum  time  of 
0.25  minutes,  a maximum  time  of  1.0  minutes,  and  a most-likely  time  of 
0.5  minutes.  This  is  the  time  expended  between  stopping  all  X and  Y 
motion  to  positioning  the  magnet  over  the  bay  and  starting  X and  Y motion 
to  transport  the  load  of  billets  to  the  work  area. 

23.  The  number  of  5-1/4  inch  billets  picked  off  of  a storage  pile  per 

unit  pick  was  entered  as  a Triangular  Distribution  having  1.0  as  the 

minimum  number  of  billets,  10.0  as  the  maximum  number  of  billets,  and 

10.0  as  the  most-likely  number  of  billets. 

24.  The  number  of  6 inch  billets  picked  off  of  a storage  pile  per  unit 
pick  was  entered  as  a Triangular  Distribution  having  1.0  as  the  minimum 
number  of  billets,  9.0  as  the  maximum  number  of  billets,  and  9.0  as  the 
most-likely  number  of  billets. 

25.  The  number  of  7-3/8  inch  billets  picked  off  of  a storage  pile  per 

unit  pick  was  entered  as  a Triangular  Distribution  having  1.0  as  the 

minimum  number  of  billets,  7.0  as  the  maximum  number  of  billets,  and 

7.0  as  the  most-likely  number  of  billets. 

26.  The  time  required  to  set  a load  of  billets  down  in  the  work  area 
begins  when  the  X and  Y travel  into  the  work  area  ends  and  continues 
until  the  X and  Y travel  away  from  the  set  down  position  begins.  This 
time  was  entered  as  a Triangular  Distribution  having  a minimum  value  of 
0.1  minutes,  a maximum  value  of  0.5  minutes  and  a mpst-likely  value  of 
0.3  minutes. 

27.  The  time  required  to  pick  up  a charge  or  a portion  of  a charge  was 
entered  as  a Triangular  Distribution  having  a minimum  value  of  0.25  minutes, 
a maximum  value  of  1.0  minutes  and  a most— likely  value  of  0.5  minutes. 

This  is  the  time  expended  between  stopping  X and  Y travel  to  positioning 
the  magnet  over  the  charge  and  the  starting  of  X and  Y travel  to  transport 
the  charge  to  the  feed  table. 

28.  The  time  required  to  set  a charge  or  a portion  of  a charge  on  the 
feed  table  begins  when  X and  Y travel  to  the  feeder  table  ends  and  con- 
tinues until  X and  Y travel  away  from  the  feeder  table  begins.  This 
time  was  entered  as  a Triangular  Distribution  having  a minimum  value  of 
0.25  minutes,  a maximum  value  of  0.75  minutes,  and  a most-likely  value 
of  0.45  minutes. 
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29.  The  time  required  to  square  up  a group  of  loose  billets  into  a 
square  array  of  tightly  packed  billets  ready  for  loading  on  the  feed  table 
was  entered  as  a Triangular  Distribution  having  a minimum  value  of  0.75 
minutes,  a maximum  value  of  2.0  minutes  and  a most-likely  value  of  1.25 
minutes.  This  is  time  expended  between  stopping  X and  Y travel  to  position- 
ing the  magnet  over  the  loose  billets  in  the  work  area  and  starting  of  X 
and  Y travel  away  from  the  work  area  to  the  next  task. 

30.  Acceleration  in  the  X direction  was  entered  as  a constant  of  1.0 
feet-per-second  squared. 

31.  Deceleration  in  the  X direction  was  entered  as  a constant  of  1.0 
feet-per-second  squared. 

32.  The  top  velocity  attainable  in  the  X direction  was  entered  as  a 
Uniform  Distribution  ranging  from  150  - 200  feet-per-minute. 

33.  Acceleration  in  the  Y direction  was  entered  as  a constant  of  1,0 
feet-per-second  squared. 

34.  Deceleration  in  the  Y direction  was  entered  as  a constant  of  1.0 
feet-per-second  squared. 

35.  The  top  velocity  attainable  in  the  Y direction  was  entered  as  a 
Uniform  Distribution  ranging  from  150  - 200  feet-per-minute. 

36.  The  size  of  a typical  5-1/4  inch  heat  was  entered  as  a constant 
of  175  billets. 


37.  The  size  of  a typical  6 inch  heat  was  entered  as  a constant  of  175 
billets. 


38.  The  size  of  a typical  7-3/8  inch  heat  was  entered  as  a constant  of 
175  billets. 

39.  The  typical  number  of  5-1/4  inch  billets  on  a railroad  car  was  entered 
as  a constant  of  60  billets. 

40.  The  typical  number  of  6 inch  billets  on  a railroad  car  was  entered 
as  a constant  of  60  billets. 

41.  The  typical  number  of  7-3/8  inch  billets  on  a railroad  car  was 
entered  as  a constant  of  60  billets. 

42.  The  time  in  minutes  a feeder  is  allowed  to  wait  before  scratching 
the  run  was  entered  as  a constant  of  20. 

Cols  4-5  (Format  12)  The  cards  required  to  accomplish  this  task 
must  be  sequentially  numbered.  The  histogram  option  is  the  only  input 
option  which  may  require  more  than  1 card  to  input  all  the  necessary 
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data  elements.  Therefore,  for  distributions  and  a constant,  the  user 
need  only  enter  a 1 in  this  field.  However,  for  histogram  data  the 
first  card  must  have  a 1 entered  in  this  field,  a 2 must  be  entered  in 
this  field  for  card  2,  and  so  on  for  as  many  cards  as  necessary. 

Cols  6-15  (Format  F10.0)  Enter  the  necessary  data  for  Field  1 as 
defined  above. 

Cols  16-25  (Format  F10.0)  Enter  the  necessary  data  for  Field  2 as 
defined  above. 

Cols  26-35  (Format  F10.0)  Enter  the  necessary  data  for  Field  3 as 
defined  above. 

Cols  36-45  (Format  F10.0)  Enter  the  necessary  data  for  Field  4 as 
defined  above. 

Cols  46-55  (Format  F10.0)  Enter  the  necessary  data  for  Field  5 as 
defined  above. 

Cols  56-65  (Format  F10.0)  Enter  the  necessary  data  for  Field  6 as 
defined  above. 

Cols  66-75  (Format  F10.0)  Enter  the  necessary  data  for  Field  7 as 
defined  above. 

Cols  76-80  Leave  Blank. 


Next  page  is  blank. 
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APPENDIX  C 


ERROR  MESSAGES 


Number  and  Description  of  Error  Messages, 

1233  - Length  of  time  between  systems  simulation  samples  is  either 

negative  or  larger  than  the  simulation  stop  time. 

1244  - Starting  time  for  taking  a sample  of  the  processing  of  the 

independent  events  is  larger  than  the  ending  time, 

1255  - Starting  time  for  taking  a system  simulation  sample  of  the 

processing  of  the  independent  events  is  negative. 

1266  - X coordinate  of  the  crane’s  starting  point  is  beyond  the 

boundaries  of  the  yard. 

1277  - Y coordinate  of  the  crane’s  starting  point  is  beyond  the 

boundaries  of  the  yard. 

1288  - Initial  value  assigned  to  the  Random  Number  Seed  is  less 

than  or  equal  to  zero, 

1299  - The  pool  level  for  one  of  the  types  of  billets  is  less 
than  the  number  of  billets  contained  in  a charge, 

1300  - The  number  of  billets  contained  in  a charge  is  a nonpositive 
number . 

1311  - The  maximum  number  of  billets  that  the  crane  can  carry 

for  one  of  the  types  of  billets  is  a nonpositive  number. 

1322  - The  X coordinate  of  one  of  the  feeders  is  beyond  the 

boundaries  of  the  yard. 

1333  - The  Y coordinate  of  one  of  the  feeders  is  beyond  the 

boundaries  of  the  yard. 

1344  - The  X coordinate  of  one  of  the  work  areas  is  beyond  the 

boundaries  of  the  yard. 

1355  - The  Y coordinate  of  one  of  the  work  areas  is  beyond  the 

boundaries  of  the  yard. 

1377  - The  initial  queue  of  heats  in  railroad  cars  waiting  to  be 

unloaded  for  one  of  the  three  types  of  billets  is  negative. 

- An  initial  value  of  one  of  the  independent  events  is 


1399 

negative. 
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1433  - The  bay  number  for  the  preceding  card  was  feasible. 

1455  - The  current  maximum  number  of  billets  allowed  in  the 
bay  card  just  read  in  is  not  a positive  number. 

1466  - One  of  the  heats  being  stored  in  the  bay  just  read  is 

either  negative  or  larger  than  the  maximum  number  of  billets  allowed. 

1488  - The  total  number  of  billets  entered  in  the  bay  just  read 

in  is  larger  than  the  maximum  number  allowed. 

1499  - The  X coordinate  of  the  bay  currently  being  read  in  is 
beyond  the  boundaries  of  the  yard. 

1500  - The  Y coordinate  of  the  bay  currently  being  read  in  is 
beyond  the  boundaries  of  the  yard. 

1544  - More  than  just  one  priority  has  been  assigned  to  the  bay 

currently  being  read  in. 

1555  - The  priority  number  assigned  to  this  bay  is  faulty. 

1577  - The  type  of  material  and  the  bay  receiving  billets  cur- 

rently being  unloaded  do  not  have  acceptable  consistent  values,  i.e,, 
they  must  both  be  zero  or  both  be  greater  than  zero  at  the  same  time. 

1599  - Either  the  type  of  material  being  unloaded,  the  bay  re- 

ceiving the  material,  or  the  heat  level  is  beyond  the  maximum  values  al- 
lowed. These  items  exceed  their  check  variables  of  MAXTYP,  MAXBAY  or 
MAXHET. 

1611  - The  bay  designated  to  receive  the  billets  queued  for  unload- 

ing does  not  accommodate  this  type  of  material. 

1644  - The  heat  level  of  the  bay  receiving  material  is  faulty. 

1655  - The  queue,  which  was  designated  to  be  unloaded,  is  empty. 

1677  - One  of  the  bays  providing  billets  to  one  of  the  feed 

tables  and  its  heat  level  does  not  have  consistent  values.  They  must 
both  be  zero  or  both  be  greater  than  zero  at  the  same  time. 

1699  - One  of  the  bays  providing  billets  to  the  feed  tables  or 

its  heat  level  exceeds  its  maximum  allowable  level.  One  or  both  of 
these  variables  exceed  their  check  variables  of  MAXBAY  or  MAXHET . 

1711  - One  of  the  bays  designated  to  providing  billets  to  a given 

feed  table  does  not  store  the  correct  type  of  billet  for  the  table  it  is 
supplying. 
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1744  - The  heat  level  of  one  of  the  bays  supplying  billets  to  a 

feed  table  is  faulty.  It  is  not  aimed  at  the  top  of  the  pile. 

1822  - A -999.0  was  not  entered  as  the  last  item  of  data  for  the 

histogram  information  previously  read  in. 

1844  - The  pointer  or  internal  reference  number  assigned  to  the 

stochastic  data  being  read  in  is  faulty. 

1866  - The  distribution  number  assigned  to  the  stochastic  data  being 

read  in  is  faulty. 

• 1899  - The  storage  area  for  stochastic  data  has  been  exceeded. 

1922  - The  card  number  is  out  of  sequence  for  the  stochastic  data 

being  read  in. 

1977  - The  Lognormal  Distribution  requested  on  the  previous  input 

card  requires  positive  non-zero  parameters* 

2000  - The  Gamma  Distribution  requested  on  the  previous  input  card 

requires  positive  non-zero  parameters. 

2022  - The  Gamma  Distribution  requested  on  the  previous  input  card 

has  a standard  deviation  which  is  too  large  or  the  mean  is  too  small 
to  yield  at  least  one  exponential  deviate. 

2044  - The  Weibull  Distribution  requested  on  the  previous  input 

card  requires  positive  non-zero  parameters. 

2077  - The  Erlang  Distribution  requested  on  the  previous  input 

card  requires  positive  non— zero  parameters  and  integer  exponential 
deviates . 

2100  - The  Chi  Square  Distribution  requested  on  the  previous 

input  card  requires  positive  non-zero  parameters  and  integer  degrees  of 
freedom. 

2133  - The  Beta  Distribution  requested  on  the  previous  input  card 

requires  positive  non-zero  parameters* 

2155  - The  Poisson  Distribution  requested  on  the  previous  input 

card  requires  positive  non-zero  parameters* 

2200  - The  Pascal  Distribution  requested  on  the  previous  input 

card  does  not  meet  one  of  the  following  requirements: 

1.  The  value  assigned  to  P must  lie  between  0 and  1. 

2.  The  value  assigned  to  K must  be  a positive  integer. 
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3.  The  minimum  observation  must  not  be  negative. 

PASCAL  - F(X)  = (K+X-1)PKQX  X = 0,  1,  2. , . AND  Q = 1-P 

( X ) 

2244.  The  Binomial  Distribution  requested  on  the  previous  input 
card  does  not  meet  one  of  the  following  requirements: 

1.  The  value  assigned  to  P must  lie  between  0 and  1, 

2.  The  value  assigned  to  N must  be  a positive  integer. 

3.  The  minimum  observation  must  not  be  negative. 

BINOMIAL  - F (X)  = (N)PXQN-X  X = 0,  1,  2 N AND  Q = 1-P 

(X) 

2288  - The  Hypergeometric  Distribution  requested  on  the  previous 

input  card  does  not  meet  one  of  the  following  requirements: 

1.  The  value  assigned  to  P must  lie  between  0 and  1. 

2.  The  value  assigned  to  N must  be  a positive  integer. 

3.  The  value  assigned  to  M must  be  a positive  integer  less 

than  N. 

4.  The  minimum  observation  must  not  be  negative. 

(NP)  (NQ  ) X *=  0,  1, NP 

HYPERGEOMETRIC  - F(X)  = (X  ) (M-X  ) M-X  - 0,  1, NQ 

(N)  Q = 1-P 

(M) 

2311  - The  distribution  requested  on  the  previous  card  requires 

a positive  standard  deviation  which  was  not  entered. 

2333  - The  distribution  requested  on  the  previous  card  has  a 

mean  value  which  lies  outside  the  minimum  - maximum  value  range. 

2355  - The  distribution  requested  on  the  previous  card  has  a zero 

or  negative  minimum  - maximum  value  range. 

2399  - The  histogram  data  just  entered  has  a probability  cell 
with  a negative  value  entered  in  it. 

2400  - The  histogram  data  just  entered  has  a probability  cell 
with  the  same  right  and  left  hand  boundaries  with  a positive  density. 

2422  - The  histogram  data  just  entered  indicates  that  a probabil- 

ity cell  has  a left  hand  boundary  which  is  larger  than  the  right  hand 
boundary . 
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2433  - The  sum  of  all  the  probability  cells  for  the  histogram  data 
just  entered  does  not  sum  to  one, 

2655  - At  least  one  of  the  bays  requires  more  heat  levels.  Expand 

check  variable  MAXHET  and  all  the  arrays  mnemonically  dimensioned  in  terms 
of  this  variable. 

2777  - A feeder  had  to  wait  for  resupplying  longer  than  the  time 

allowed  as  per  stochastic  input  number  42. 

2811  - At  least  one  of  the  heat  queues  requires  more  heat  levels. 

Expand  check  variable  MAXQUE  and  all  the  arrays  mnemonically  dimensioned 
in  terms  of  this  variable. 

2955  - The  yard  is  cleaned  out  of  one  of  the  three  different 

types  of  billets.  The  error  number  listed  after  this  one  indicates  the 
type  of  billet  (1  = 5-1/4,  2=6,  3 = 7-3/8)  the  yard  is  out  of. 


Next  page  is  blank. 
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APPENDIX  D 


COMPUTER  LISTING  OF  PROGRAM 
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0064  IF {X.GT.SMAX (5) ) SMAX(5)  = X GLM10820 

0065  CALL  SNAP  GLM10830 

0066  RETURN  GLM10840 

0067  END  GLM10850 
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APPENDIX  E 


LETTER  OF  REQUEST  FOR  THE  STUDY 
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DEPARTMENT  OF  THE  ARMY 
SCRANTON  ARMY  AMMUNITION  PLANT 
SCRANTON.  PENNSYLVANIA  18501 


IN  REPLY  REFER  TO: 

SARSC  - c 


19  December  1975 


SUBJECT:  Billet  Yard  Crane 


Commander 

US  Army  Armament  Command 

ATTN : ANSAR'  S A 

Rock  Island,  IL  61201 


1.  Reference  is  made  to  the  visit  of  Mr.  Jerry  Moeller  to  Scranton 
AAP  relative  to  preparing  a mathematical  model  of  the  Billet  Yard 
Crane  operation, 

2.  Subsequent  to  reference  visit  and  further  discussions  thereof, 
this  letter  formalizes  our  request  for  support. 

3.  It  is  understood  that  a mathematical  model  can  be  structured  but 
that  additional  operational  data  is  required  for  satisfactory  input. 
Your  office  is  requested  to  provide  the  model.  When  adequate  data  is 
collected,  a further  review  of  action  will  be  made. 

4.  Your  assistance  in  this  matter  is  greatly  appreciated. 
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DISTRIBUTION  LIST 


No.  of  Copies 


11 

2 

1 

Commander 

US  Army  Armament  Command 
ATTN:  DRSAR-SAL 

DRSAR-SAS 
DRSAR-PPI-WM 
Rock  Island,  XL  61201 

1 

Commander 

Rock  Island  Arsenal 
ATTN:  SARRI-LPL 

Rock  Island,  IL  61201 

1 

Commander 

Badger  Army  Ammunition  Plant 
Baraboo,  WI  53913 

1 

Commander 

Cornhusker  Army  Ammunition  Plant 
Grand  Island,  NE  68801 

1 

Commander 

Gateway  Army  Ammunition  Plant 
St.  Louis,  MO  63143 

1 

Commander 

Hays  Army  Ammunition  Plant 
Pittsburgh  (Hays) , PA  15207 

1 

Commander 

Hols ton  Army  Ammunition  Plant 
Kingsport,  TN  37662 

1 

Commander 

Indiana  Army  Ammunition  Plant 
Charleston,  IN  47111 

1 

Commander 

Iowa  Army  Ammunition  Plant 
Burlington,  IA  52601 

1 

Commander 

Joliet  Army  Ammunition  Plant 
Joliet,  IL  60436 

1 

Commander 

Kansas  Army  Ammunition  Plant 
Parsons,  KS  67357 
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DISTRIBUTION  LIST  (cont'd) 


Commander 

Lake  City  Army  Ammunition  Plant 
Independence,  MO  64050 

Commander 

Lone  Star  Army  Ammunition  Plant 
Texarkana,  TX  75501 

Commander 

Longhorn  Army  Ammunition  Plant 
Marshall,  TX  75670 

Commander 

Louisiana  Army  Ammunition  Plant 
Shreveport,  LA  71102 

Commander 

Milan  Army  Ammunition  Plant 
Milan,  TN  38358 

Commander 

Newport  Army  Ammunition  Plant 
Newport,  IN  47966 

Commander 

Radford  Army  Ammunition  Plant 
Radford,  VA  24141 

Commander 

Ravenna  Army  Ammunition  Plant 
Ravenna,  OH  44266 

Commander 

Riverbank  Army  Ammunition  Plant 
Riverbank,  CA  95367 

Commander 

St.  Louis  Army  Ammunition  Plant 
St.  Louis,  M0  63120 

Commander 

Scranton  Army  Ammunition  Plant 
Scranton,  PA  18502 


DISTRIBUTION  LIST  (cont'd) 


No.  of  Copies 

1 Commander 

Sunflower  Army  Ammunition  Plant 
Lawrence,  KS  66044 

1 Commander 

Twin  City  Army  Ammunition  Plant 
Minneapolis,  MN  55112 

1 Commander 

Volunteer  Army  Ammunition  Plant 
Chattanooga,  TN  37401 

12  Defense  Documentation  Center 

Cameron  Station 
Alexandria,  VA  22314 
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